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INTRODUCTION 


Estimation of solid state electrostatic potentials, fields and field 
gradients In Ionic solids has been of considerable Interest during the 
past few years. Several methods have been developed starting from a simple 
point charge model of the solid to an elaborate multi pole-model of the 
lattice^”®. The results are of varying degrees of success depending on the 
accuracy of the a) crystal structure data b) dipole and quadrupole 
polarizability values and c) the number of Ions of the lattice taken Into 
consideration In arriving at the potential, though there exists considerable 
controversy regarding the validity of the multipole model of the lattice, 
and the usefulness of such calculations due to the uncertainties In the para- 
meters (a) and (b) mentioned above, and the appropriateness of allocating 

a certain charge value to the lattice site, and lastly the uncertainty In the 

T 8 

Sternhelmer antlshlelding factor * , the method has yielded satisfactory 

9-12 

results In a number of cases 

The main purpose of our work is to develop a straightforward self- 
consistent method, that can easily be understood, without presenting 
serious mathematical conceptual problems, by the senior undergraduate 
students who are familiar with basic principles of electrostatics. 

The method will be a direct practical application of basic electrostatics 
to solid state and also helps in the understanding of the principles of 
crystal structure. This report presents the necessary mathematical 
equations derived from first principles and the systematic computational 
procedure developed to arrive at the solid state EF 6 (electrostatic 
field gradients) values. No such calculations have so far been made for 
the oxides of third and fifth group elements except for Al202^ which is 
used for testing our method and computer program. The program will be 
applied to estimate the EFG values (or ^ values) in the oxides of 
gallium, iridium, arsenic, antimony and bismuth. 
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Computational procedure 
This involves;- 

a) Environmental search. determination of the cartesian co-ordinates of 
all atoms within a chosen distance from the ion at the lattice 
site of interest. 

b) Evaluation of the solid state electrostatic potential (V), field 

(Vj =1^.1 = X. y. z) and field gradient (Vi,. = a^V , i. j = x, y, 2 ) 
at the lattice sites of interest. 

c) Self consistent evaluation of the induced dipole and quadrupole 
moments of all polarisable ions of the lattice. 

d) Final evaluation of the EFG at the site of interest due to the charge, 
dipole and quadrupole nature of all the lattice ions. 

environmental search 

starting with the crystal structure data consisting of the unit cell 
dimensions A^, B^, C^ the number of atoms per unit cell H).. and their 
position co-ordinate parameters U. V. W. (Input data for the program), 

co-ordinates of all atoms within x(= 1,2,3 ) unit cells aVfay 

in the three dimensions are located in the crystalline axis system. 

The progress of the search sequence is shown in Fig. 1. Starting with 
the given crystalline co-ordinates A^U, B^V, C^W of any one of the N^ atoms 
in the unit cell, all similar atoms in neighbouring unit cells are located 
as Indicated in figure 1 for x = 1. There will be 27 identical atoms if 
the search is confined to one unit cell aVfay, making up the total number 
of ions included in the computation to 27 N^. If the search is extended to 

3 

X unit cells the total number of ions will be N^ * (2x + 1) N^. The choice 
of the sphere of influence being in multiples of unit cells is considered as 
more appropriate than an arbitrary choice of R value as 50A°U or 100 A°U, 
because charge neutrality is preserved in units of crystalline unit cells. 

It will be necessary to convert the crystalline co-ordinates into cartesian 
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co-ordinates if the unit cell symnetry Is hexagonal or rhcMtbchedral 
using the following transformation matrices. 


Hexagonal to cartesian : 


Rhombohedral to cartesian : 


sine sina/. 


sine - si no 


/2 



cose 




cose 


- sine 


cose 


where a is the rhombohedral angle and e is the angle between z-axis 
(axis of symmetry) and the crystalline axes. Sine = 

The first part of the program computes the cartesian co-ordinates of all 
the atoms relative to an arbitrary origin and stores them in the X^, 
^0* ^0 origin is shifted to one of theiipnequi valent 

sites of interest and the relative co-ordinates and distances of all the 
ions are computed and stored in another working array XR, YR, ZR, RR for 
computing the potential, field, and field gradient at that site. 

Potential, Field and Field Gradient 

The expressions for fields (V. = |j , (i = x, y, z) and field 
gradients , i , j , = x, y, z) at a point P (x, y, z) in 




space due to a multipole situated at the origin having a point charge C , 
dipole moment D (= i D + j D + k D_) and quadrupole moment 

A y z 


^ -- 


arp obtained by direct differentiation of the potential expressions for 
the monopole (VM), dipole (V^), and quadrupole (V^) * 



since the general expressions derived In cartesian co-ordinates could 
not be reduced Into a very compact form they are presented In an analytical 
form Indicating the Individual contributions of each dipole component and 
quadrupole tensor element to the field and field gradient. 
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The total potential (V) and hence the field (V^\ and field gradient 

(V 4 .) at any Ion site s Is . ^ 

^ Kit Pi « 

V(sj = V, -V 


where n Includes all of the Ions of the lattice. But they cannot be 
computed without the knowledge of the multlpole moments which 

are, in turn, dependent on and 
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of the lattice as given by the equations above. As such, a self- 


consistent solutibi. Is sought for the estimation of the E.F.G's. As can be 
seen from the detailed expressions, the coefficients of the multipole moments 
of the lattice charge are entirely determined by their position co-ordInateS 
only. Knowing the crystal structure data, the environmental search part 
of tne program computes the cartesian co-ordinates of all the lattice charges 
relative to the chosed Ion site and stores them In the XR, YR, ZR, RR array. 
Hence the coefficients can be computed Independently. 

For computational convenience the detailed expressions given above are 
expressed in terms of a set of and P 4 functions. By using these functions 
the development of the computer program has become considerably simpler since 
r.iany of the coefficients repeat themselves, “Jo save the computational time on 
the computer the coefficients are presented in the form of a matrix (C^j). 

The number of computations needed for each Ion considered are reduced by 50% 
by this method. The S and P functions and the C matrix in terms of S and P 


functions |s given below. Only the elements given by the S functions need 
to be computed. If an element (CIO) occurs In a subsequent coluimi as Ckl 
(CKL element) It Is given as CIO. 
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Self consistent solutions of muUipole moments 

Taking one by one, each of the atoms chosen to be Included In the 
computation the 90 elements are computed for each atom and the final sum 
Is stored In the term matrix 


T.j= 2 c, 

ns I 


The first three rows of the T-roatrIx represent the contribution of the 
lattice to the electrostatic field at the Ion s1te^)vfhen summed along the 
row after multiplication by the corresponding muUipole moment D. (1 » 1, tO 
monopole moment or charge Cp; Dg, Dj, D^, components of the 

Induced dipole moment; Dg, Dg, D^, Dg, Dg. D^g » 

Induced quadrupole moment tensor elements respectively) and the last six 
rows the electrostatic field gradient. Thus the rows of the C/T matrix give 


t a 1 


V ■ 
Xxv 

i * 2. 

Vyy i 

Vyt 
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To start with an arbitrary set of D-values (charge, dipole and 
quadrupole moments) Is assigned and provided as Input data to the computer 
along with the values of dipole and quadrupole polarlsabi 11 tleso^ and^. 
Follov.'ing the Gauss-ZIllard iteration process the self consistent Dipole 
and Quadrupole moments are obtained. Since there are only 9 independent 
equations, three for t.ie fields Vy, and six for the field gradients 
Vxx* Vyy, V,^. Vy^, 006 COO 60^6 for D^,Dy, and Qyy, 

Qyy* ^zx* charge value (Cp-D^) should be assigned b.ased on the 

ionic nature of the solid or one can determine the best possible value so as 
to be consistent with the experimental value of field gradients when they are 
known. The 9 equations are represented by two sets 
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Note D| Is eliminated from the equations by incrementing the index 
on the left by one. Inserting the assigned values for D^.'s in the first 
equation (i - 1) a new value is obtained for 02* i.e. D^. This value 
of D„ is used in solving the second equation (i = 2) for D, i.e. D„. 

A ^ J ' 

In the third equation (i « 3) new values of and D are used and D, 
i.e. 1^ is obtained and the process continues using improved values 
progressively until the 9-D values are obtained. After completing the first 
cycle the new D-values are compared with the old values and the cycle is 
repeated using the new values until the total difference between successive 
iterations is ^0.0001. 
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In a typical case this accuracy was achieved in eight iterations. 

If the lattice consists of more than one inequivalent ions, either due to 
their chemical or physical nature the number of unknowns becomes 9S, where S 
is the number of inequivalent lattice sites where polarisable charge 
distribution exists. In general metallic atoms are treated merely as point 
charges, since their polarisability values are negligible compared with 
non-metal lie atoms. The T-matrix should be evaluated at each of the S-sites 
to render the 9S simultaneous equations. 




EVALUATION OF THE EFG AT THE QUADRUPOLE NUCLEUS 


The final phase of the prog»"am consists almost of a repetition of 
the above computation up to the evaluation of the T-matrix, starting with 
the lattice site of Interest* v^here the quadrupole nucleus Is located* 
as origin. Using the original array XO* YO* ZO of the position co*ord1nates, 
the cartesian co-ordinates of all atoms relative to the quadrupole 
nucleus of interest are computed:- 

XR * X0(Q) - X0(n) etc. and so on and the program proceeds until the term 
matrix T^j is obtained. Using the self consistent multi pole moments 
evaluated in the previous phase* the electrostatic fields and field gradients 
are readily obtained from 
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1=1-3 Fields and i = 4 - 9 Field gradients. 

The final output is obtained in the form of 9 x 13 matrix. The 

1st column indicates the contribution due to the lattice charges only V 

dipole 

2nd, 3rd and 4th due to individual /components Vpy* Vp^* 5th 

column total contribution due to the dipoles Vp* 6 - 11th individual 

contributions due to the six quadrupole tensor elements 

^0 * ^0 ‘ '^0 ’ ’ '^0 » column the total contribution 

^xx ^yy ^zz ^xy ^z ^zx 

of quadrupoles Vq and finally the 13th column is the total contribution 
due to all multi poles Vj. The program is tested by computing the EFG 
values for AI 2 O 2 for which data is available^, and is presently being 
extended to study the oxides of arsenic* antimony, bismuth, indium and 
gallium. 

Experimental nuclear quadrupole resonance data is available for the V 
group oxides (As, Sb and Bi). Search for NQR signals in the oxides of iKldluro 
and gallium has so far not yielded any results at room temperature. These will 
be re-investigated at liquid nitrogen temperatures at the Queen Elizabeth 
College NQR laboratories. 
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